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Spin Gap and Superconductivity in Weakly Coupled Ladders: 
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Effects of the interladder one-particle hopping, t±, on the low-energy asymptotics of a weakly 
coupled Hubbard ladder system have been studied, based on the perturbative renormalization- 
group approach. We found that for finite intraladder Ifubbard repulsion, U, there exists a 
crossover value of the interladder one-particle hopping, t±c- For < < txc, the spin gap 
metal (SGM) phase of the isolated ladder transits at a finite transition temperature, Tc, to the 
d-wave superconducting (SCd) phase via a two-particle crossover. In the temperature region, 
T < Tc, interladder coherent Josephson tunneling of the Cooper pairs occurs, while the inter- 
ladder coherent one-particle process is strongly suppressed. For t±c < t±, around a crossover 
temperature, Tcrosa, the system crosses over to the two-dimensional (2D) phase via a one-particle 
crossover. In the temperature region, T < Tcross, the interladder co/iereni band motion occurs. 

KEYWORDS: doped Hubbard ladders, interladder coupling, one-particle crossover, two-particle crossover, d-wave 
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Last year Uehara et al.n discovered the superconduc- 
tivity signal with Tc — 12K in the doped spin ladder, 
Sro.4Ca13.6Cu24O41.84, under a pressure of 3GPa. The 
electric properties of the compound arcL determined by 
the hole-doped ladders instead of chains .0^ A remarkable 
feature of the. doped ladder is the existence of a spin ex- 
citation gap.cP The superconducting transition under a 
high pressure suggests that interladder one-particle hop- 
ping induced by the applied pressure plays an important 
role. Recent experiments on the resistivity alonaihe lad- 
der, pc of the single crystal Sr2.5Caii.5Cu2404iQ' shows 
that the superconductivity sets in below 10 K under 
3.5GPa ~ 8GPa with the temperature dependence of pc 
changing gradually from T-linear to T^. The anisotropy 
of the resistivity also indicates the dimensional crossover 
from ID to 2D with increasing an applied pressure. 

In this paper, to elucidate the nature of superconduc- 
tivity in the doped ladder^compound under pressure, we 
consider Hubbard laddersB coupled via a weak interlad- 
der one-particle hopping. In the case of the isolated Hub- 
bard ladder, the most relevant phase is characterized by 
a strong coupling fixed point and is denoted by "phase II' 
by FabriziocP and "CISO phase" by Balents and Fisher. EP 
In tJpis jihase, only the total charge mode remains gap- 
IcsgS'El'Erlj'' and consequently, the c?-wave superconduct- 
ing correlation becomes the most doniinanLonc as long 
as the intraladder correlation is weakfl'Lf B'El-' From now 
on we call this strong coupling phase "spin gap metal 
(SGM) phase" . So far, the effects of interladder hopping 
on the SGM phasejiave been studied through mean fiekl 
approximationstj'Ej' and power counting arguments. IliP 

The central problem here is how a weak interladder 
one-particle hopping, t±, affects the low-energy asymp- 
totics of the system. Based on the perturbative renormali- 

* E-mail:kishine@ims. ac.jp 



zation-group (PRG) approach, we here study one-particle 
and two-particle crossovers JIS' when we switch on t± and 
the intraladder Hubbard repulsion, U, as perturbations 
to the system, specified by the intraladder longitudinal 
(transverse) hopping, t{t'): U,t± ^ t,t'. A similar ap- 
proach has been considered for^the problem of a cou- 
pled chain system by Boies et a/iilP The intraladder one- 
particle process is diagonalized in terms of the bonding 
(B) and antibonding (A) bands. As shown in Fig. 1(a), 
we linearize the dispersions along the legs on the bond- 
ing and antibonding Fermi points, ±fci?m(m = A, B). In 
Fig. 1, R and L denote the right-going and left-going 
branches, respectively. The Fermi velocities in principle 
depend on the band index as vpm — 2t smkpj^ but we 
assume throughout this work that vprn = vp and drop 
the band index, since the difference in the Fermi veloc- 
ities does not affect the asymptotic nature of the SGM 
phase at least for small i'/t.BQ) In all the four branches 
(LB, LA^ RA, RB) of linearized bands, the energy vari- 
ables, e^m {v — R,L ; m ~ A, B) run over the region, 
—E/2 < e^m < E/2, with E denoting the bandwidth 
cutoff. 

The intraladder Hubbard repulsion generates the scat- 
tering processes depicted in Fig. 1(b). The processes 

are specified by dimensionless coupling constants gi}^ 

(2) . . 

and gli denoting backward and forward scatterings, re- 
spectively, with the flavor indiceaj-* p — 0, f,t denoting 
intraband scattering, interband forward scattering and 
interband tunneling processes, respectively. The usual 
coupling constants with a dimension of the interaction 
energy are . We neglect the interband back- 

ward processes such as {RB RA; LA LB) on the 
grounds that these processes do not seniiiusly modify the 
asymptotic nature of the SGM phase. ErLl^ 

When we switch on gfi and tj^ as perturbations to 
the system with t and t\ and the temperature scale 
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Fig. 1. (a) Four hvanches(LB, LA, RA, RB) of linearized bands 
with the bandwidth cutoff E and (b) intraladder two-particle 
scattering vertices 3^'' . The solid and broken lines represent the 
propagators for the right-moving and left-moving electrons, m 
and m denote different bands. 



Fig. 2. Schematic illustrations of the one-particle and the two- 
particle process (in the case of d-wave superconductivity chan- 
nel). In the one-particle process, a particle hops from one ladder 
to a neighboring one, while in the two-particle process, a pair of 
particles hops from one ladder to a neighboring one. 



decreases, two kinds of dimensional crossovers- a^one- 
particle crossover and a two-particle crossover ,EJ'E3' oc- 
cur. We illustrate one- and two-particle processes in 
Fig. 2. In the case where the interladder one-particle 
hopping modified by the intraladder self-energy becomes 
the most relevant, a particle hops coherently from one 
ladder to a neighboring one and the system crosses over 
to a two-dimensional system via the one-particle crossover. 
In the case where the interladder two-particle process 
generated by the mierladder one-particle hopping and in- 
iraladder two-particle scattering processes becomes the 
most relevant, a pair of composite particles hops coher- 
ently from one ladder to a neighboring one and the sys- 
tem crosses over to a long-range-ordered phase via the 
two-particle crossover. In Fig. 2, we show only the two- 
particle hopping in the d-wave superconducting channel, 
which corresponds to the interladder Josephson tunnel- 
ing of Cooper pairs. 

To study the competition between the one-particle crossov|;c 
and the two-particle crossover, we set up scaling equa- ' 
tions for the interladder one-particle and two-particle 
hopping amplitudes and study their low-energy asymp- 
totic behavior. We start with the path integral repre- 
sentation of the partition function of the system, Z = 

Ve^ , where the action consists of four parts, 



tion will be renormalized as I goes from zero to infinity. 
The scaling equations for all the processes studied here 
are depicted in Fig. 3. Diagrammatic expansions for 



3. Diagrammatic representations of the scaling equations for 
the intraladder scattering vertices(a), intraladder one-particle 
propagator (b), interladder one-particle hopping amplitude(d), 
and interladder two-particle hopping amplitude(d). A black cir- 
cle, zigzag line, and shaded square represent the intraladder 
two-particle scattering processes (a combination of g^'^), the 
interladder one-particle hopping amplitude, t±, and the in- 
terladder two-particle hopping amplitude, ySCd respectively. 
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where 5*11 , 5*11 , and denote the action for 

the intraladder one-particle hopping, intraladder two- 
particle scatterings, interladder one-particle and inter- 
ladder two-particle hopping, respectively. V symbol- 
izes the measure of the path integral over the fermionic 
Grassmann variables. 

The idea of scaling is to eliminate the short-wavelength 
degrees of freedom to relate effective actions at successive 
energy scales. Based on the bandwidth cutoff regulariza- 
tion scheme, we parametrize the cutoff as E{1) = -Boe ' 
with the scaling parameter, I, and study how the ac- 



the intraladder scattering processes up to the 3rd or- 
der are given in Fig. 3 (a). After taking account of the 
field rescaling procedure which originated from llie scal- 
ing of the intraladder one-particle propagatorSjEj-* given 
in Fig. 3(b), we obtain the appropriate scaling equations 
for gj^^ . Full expressions of the scaling equations for g};*' 
are found by setting g^^^ — gf^ = in Eq.(A5) of rcf. [6]. 
Starting with the Hubbard-type initial condition 



9« 



(0) = [7 EE U/Attvf > 0, 



(2) 
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the scaling equations give the fixed point 



9i''* - -1, 

(2)* _ 3^20 
yo — 4 : 



9?* = 0, 

(2)* _ 1+2U 
9f - 4 : 



(1)* _ 1 

V" ' (3) 
9t 



1. 



Henceforth, we take U = 0.3 for illustration. For any 
U > 0, the results are similar to those shown below. We 
show the scaling flows for g^*"* in Fig. 4(a) as functions of 
the scaling parameter, /. These flop^ scale the isolated 



Hubbard ladder to the SGM phase.l 

The one-particle crossover is determined through the 
scaling equation 

- 1 - (5^^ + 9^ + 9?' + 9?' + 9i''' 

-^<^'-<^<^-</f<^~<^<^l (4) 

which is represented by Fig. 3(c). It is seen from (4) 
together with the fixed point (3) that 



d\iitj_{l)/dl 



- 7/8, 



(5) 



and consequently t±{l) becomes always irrelevant at the 
final stage of the scaling procedure. However, t±(l) grows 
at an early stage of scaling since the intraladder cou- 
plings do not grow sufficiently as yet (see Fig. 4). Once 
i!;j_(/) = t±{l)/Eo attains an order of unity, the weakly 
coupled ladder picture breaks down and the system is 
scaled to aJj two-dimensional" system via the one-particle 
crossover.Ej^ In the upper planes of Figs. 4(b) and 4(c), 
we show the scaling flows of t±{l) with the initial con- 
ditions, ij_{0) = ~ 0.01 and 0.04, respectively. We 



Fig. 4. Scaling flows of g^,'^(a), t^(Z)(the upper half planes of 
(b),(c)), and ^^'-''^(the lower half planes of (b),(c)). For the re- 
gion, I > Min(/c, 'cross), the scaling flows, denoted by broken 
curves, have no physical meaning, since the weak coupling pic- 
ture breaks down in the region. 



see that for i±o = 0.01, i±{l) never reaches unity and 
one-dimensional crossover never takes place, while for 
i±o = 0.04, i±{l) exceeds unity at the scaling parameter 
'cross which is defined by 



1. 



(6) 



Thus we specify the one-dimensional crossover by Ic 
although it merely has qualitative meaning. 



Interladder two-particle processes are decomposed into 
CDW, SDW, SS (singlet superconductivity) and TS (triplet 
superccmductivity) channels as in the case of the coupled 
chains.113' In this case, there are additional flavor indices, 
fj, — 0, f,t for each channel. Then the two-particle hop- 
ping amplitudes are specified as V^f^ (Af=CDW, SDW, 
SS, TS; n = 0,f,t). The SS channel can be decom- 
posed into the s-wave spuerconductivitv. [iSCf ) and d- 
wave superconductivity (SCd) channels .&Q'&B' For the 
SCd channel, the action for the interladder two-particle 
hopping is written in the form 



S 



(2) 

_LSCd 



^)^F^^^O*cd(Q)Oscd(Q). (7) 



Q 



The SCd pair-field is given by OscdiQ) = 0^^{Q) - 
0^s^(g), where O^^^' (Q) = P'^/'Ek,.^ RmA-K + 
Q)Lm'.-a-{K), with Rm,cr{Lm,a) being the Grassmann 
variable representing the right (left)-moving electron in 
the band m with spin a and K — (fc||, fcj^,i£„) and Q — 
{q\\, qj_, itj„) with the momentum along the leg and rung 
being specified by || and _L, respectively, and fermion 
and boson thermal frequencies, e„ — {2n + 1)t:/ (3 and 
ujfi — 2mr/f3, respectively. 

The lowest-order scaling equation for V^'~^'^ is depicted 
in Fig. 3(d), and is written as 



dl 



^-[h{i)9'^''{i)y 



1 



+ 25SCd(;)^scd(;)_^ [^scd(;)]^^ 



where the transverse momentum transfer of the pair is 
set a.t q± — 0. The coupling for the SCd pair field is given 

by 5'°' = ^(^?^^ + 5f ^ - ff^'^ - 5^'^ By lengthy but 
straightforward manipulation, we obtain similar scaling 
equations for all V.f'{l). We have solved them with the 



initial conditions 



(9) 



and confirmed that y^^^d g^i-^^g^yg dominates the other 
channels. This situation is quite reasonable on the physi- 
cal grounds that the interladder pair tunneling stabilizes 
the most dominant intraladder correlation, i.e, d-wave 
superconducting correlation. Below, we focus on the d- 
wave superconducting channel. The third term of the 
r.h.s of eq.(8) causes divergence of ]/^c:d critical 
scaling parameter Ic determined by 



V 



SCd 



(Ic) = -00 



(10) 



In the lower half planes of Figs. 4(b) and 4(c), we show 
the scaling flows of V^^'^H) for i^Q = 0.01 and 0.04. 

Figures 4(b) and 4(c) show that, in the case of i±Q = 
0.01, only the two-particle crossover occurs, while in the 
case of t±o = 0.04, the one-particle crossover dominates 
the two-particle crossover. We identify the scaling pa- 
rameter with the absolute temperature as / = In . 
Thus, based on eqs.(6) and (10), we define the one- 
particle crossover temperature, Tdoss, and the d-wave su- 
perconducting transition temperature, Tc, as 



cross — ^ 



(11) 
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Tcross and Tc correspond to Txi and li^, respectively, in- 
troduced by Bourbonnais and CarontS' for weakly cou- 
pled chains. 

By solving the scaling equations for various tj^o with 
the fixed value of the Hubbard repulsion, U = 0.3, we 
obtained the phase diagram of the system in terms of 
t^o and the reduced temperature T — T/Eq, as shown 
in Fig. 5. Roughly speaking, we can regard increasing t±Q 



cross to the SGM 
cross, the inter- 



Fig. 5. Phase diagram of tiic weakly coupled Hubbard ladder sys- 
tem. SGM, SCd and 2D denote the sping gap metal phase, the 
d-wave superconducting phase and the two-dimensional phase, 
respectively. The interladder one-particle hopping and temper- 
ature are scaled by the initial bandwidth cutoff, Bo- 



as applying the pressure under which the bulk supercon- 
ductivity was actually observed. We found that there 
exists a crossover value of the interladder one-particle 
hopping, i±c ^ 0.025. 

For < t_Lo < t_Lc, the phase transition into the d-wave 
superconducting (SCd) phase occurs at a finite transition 
temperature, Tc, via the two-particle crossover. In the 
temperature region, T < Tc, coherent Josephson tun- 
neling of the Cooper pairs in the interladder transverse 
direction occurs. Here we have to be careful about iden- 
tification of the finite temperature phase above Tc, where 
the system is in the isolated ladder regime. As the tem- 
perature scale decreases, the isolated ladders are grad- 
ually scaled toward their low-energy asymptotics, the 
SGM phase. The gradual change of the density of dark- 
ness in the SGM phase in Fig. 5 schematically depicts this 
situation. The SGM phase is characterized by the strong 
coupling values of intraladder couplings, g^ = gp'' = 1 
and ^Q^^ = — 1, (see (3)). The critical scaling parameter, 
Ic, is in the region, Ic > 5.3, around which the intral- 
adder coupling constants almost reach thier fixed point 
values (see Fig.4.(a)). Thus we expect that the spin gap 
is well developed near Tc. Within the framework of the 
PRC approach, we cannot say for certain whether the 
spin gap survives in the SCd phase or not. 

For t±c < i_LOi around the crossover temperature Tcross, 
the system crosses over to the 2D phase via the one- 
particle crossover. The crossover value of the scaling 
parameter, /cross, is in the region /cross < 4.6, around 
which the intraladder coupling constants are far away 
from thier fixed point values (see Fig.4(a)). Thus it is 



disputable to assign the phase above Tc 
phase. In the temperature region, T < T^ 
ladder coherent band motion occurs. Then the physical 
properties of the system would strongly depend on the 
shape of the 2D Fermi surface. 

Finally we briefiy compare the present case-with the 
case of coupled chains within the PRG scheme .113' The 
scaling equation for the interchain one-particle ho 
instead of (4), gives Tcross ~ [t_Lo/£^o]^^^^~^' 
where one has the exact result for the anomalous ex- 
iupnent, 6 < 1/8, for the non-half-filled Hubbard model. 
O Consequently tj_ becomes always relevant, contrary 
to the coupled ladder case. This difference reflects the 
fact that the isolated Hubbard ladder belongs to the 
strong couphng universality class (SGM phase), while 
the isolated Hubbard chain belongs to the weak coupling 
(Tomonaga-Luttinger) universality class. In that sense, 
we can conclude that the spin gap opening strongly sup- 
presses the one-particle crossover so that the d-wave su- 
perconducting transition via the two-particle crossover is 
strongly assisted in the weakly coupled ladder system. 
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